Determination of scattering lengths from measurement of n^n atom lifetime 
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Abstract 



The DIRAC experiment at CERN has achieved a sizeable production of n^n atoms and has significantly improved the precision on 
■ ^ its lifetime determination. From a sample of 21227 atomic pairs, a 4% measurement of the S-wave nn scattering length diff'erence 

^ |ao - fl2l = (0-25331[j ™^°[j^j -ooovsLyst) attained, providing an important test of Chiral Perfiirbation Theory. 
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1. Introduction 

Pionium {A2jt) is the tt^tt" hydrogen-like atom, with 378 fm 
Bohr radius, which decays predominantly into tt'V | 1 |. The al- 
ternative yy decay accounts for only ~ 0.4% of the total rate ||2l. 
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Its ground-state Ufetime is governed by the nn S-wave scatter- 
ing lengths fl/, with total isospin / = 0, 2 |[Tl[3l: 

^2n0 = \ P*(flO - fl2)'(l + ml., (1) 

where p* - ^M^+ - - (l/4)a^Ml+ is the jr" momentum 
in the atom rest frame, a is the fine-structure constant, and 6 - 
(5.8+1 .2) ■ 10 ~ is a correction of order a due to QED and QCD 
O which ensures a 1% accuracy of equation ([T]l. The value of 
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do and a2 can be rigorously calculated in Chiral Perturbation 
Theory (ChPT) \A. 5 1, predicting a,) - ^2 = (0.265 + 0.004)M^,i 
and the Hfetime r = (2.9 ± 0.1) ■ 10"'^ s |6|. The mea- 
surement of T2,fi provides an important test of the theory since 
flo - (32 is sensitive to the quark condensate defining the spon- 
taneous chiral symmetry breaking in QCD |7|. The method re- 
ported in this article implies observation of the pionium state 
through its ionization into two pions. Given its large Bohr 
radius, this is directly sensitive to kti scattering at threshold, 
Mj,„ ~ 2M„+, and thus delivers a precision test of the theory 
without requiring threshold extrapolation, as for semileptonic 
Ke4 decays [8l, or substantial theoretical input as for — > Stt 
decays |i9j. 

2. Pionium formation and decay 

In collisions with target nuclei, protons can produce pairs 
of oppositely charged pions. Final-state Coulomb interaction 
leads to an enhancement of n^n^ pairs at low relative cm. mo- 
mentum (Q) and to the formation of A2„ bound states or pio- 
nium. These atoms may either directly decay, or evolve by ex- 
citation (de-excitation) to different quantum states. They would 
finally decay or be broken up (be ionized) by the electric field of 
the target atoms. In the case of decay, the most probable chan- 
nel is iTT^TT*' and the next channel is 77 with a small branching ra- 
tio of 0.36%. In the case of breakup, characteristic atomic pion 
pairs emerge |10|. These have a very low Q (< 3 MeV/c) and 
very small opening angle in the laboratory frame (< 3 mrad). 
A high-resolution magnetic spectrometer (Ap/p ~ 3 ■ 10"^) is 
used ifTTI to split the pairs and measure their relative momen- 
tum with sufficient precision to detect the pionium signal. This 
signal lays above a continuum background from free (unbound) 
Coulomb pairs produced in semi-inclusive proton-nucleus in- 
teractions. Other background sources are non-Coulomb pairs 
where one or both pions originate from a long-lived source 
(j], 77', A, . . .) and accidental coincidences from difl'erent proton- 
nucleus interactions. 

The first observation of A2„ was performed in the early 
1990s IT2I. Later, the DIRAC experiment at CERN was able 
to produce and detect ~ 6000 atomic pairs and perform a first 
measurement of the pionium lifetime 1131 . We now present fi- 
nal results from the analysis of ~ 1.5-10^ events recorded from 
2001 to 2003. Compared to the results in 1 13 J, this analysis has 
reduced systematic errors and improved track reconstruction, 
mostly due to the use of the GEM-MSGC detector fTP\ infor- 
mation, which leads to a larger signal yield. The present data 
come from collisions of 20 and 24 GeV/c protons, delivered by 
the CERN PS, impinging on a thin Ni target foil of 94 or 98 yum 
thickness for different run periods. 

3. Pionium detection and signal analysis 

Low relative-momentum prompt and accidental ;r^7r" pairs 
are produced at the target and selected by the multi-level trig- 
ger when their time difference, recorded by the two spectrom- 
eter arms, is |Af| < 30 ns. A suitable choice of the target ma- 
terial and thickness provides the appropriate balance between 



the A2n breakup and annihilation yields, with reduced multiple- 
scattering |[T4][15]. For a thin Ni target, of order ~ 10"^^ Xq, 
the relative cm. momentum Q of the atomic pairs is less than 
~ 3 MeV/c and their number is ~ 10% of the total number of 
free pairs in the same Q region lfT6l . The experiment is thus 
designed for maximal signal sensitivity in a very reduced re- 
gion of the tt^t:^ phase space. This is done by selective trigger- 
ing and by exploiting the high resolution of the spectrometer 
and background rejection capabilities. The longitudinal (Ql) 
and transverse (Qt) components of Q, defined with respect to 
the direction of the total laboratory momentum of the pair, are 
measured with precisions 0.55 MeV/c and 0.10 MeV/c, respec- 
tively. 

The double differential spectrum of prompt tt^tt" pairs A^pr 
(defined by |Af| < 0.5 ns), composed of atomic ha, Coulomb 
A^c, non-Coulomb A^nC, and accidental A^a^c pairs, can be x^- 
analysed in the (Qt, Ql) plane by minimizing the expression 
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[M'J - F'^ - F-^] 
[M'j + (cr'l)^ + (tr^)2] 



Here 



M(Qt, Ql) = 



AQtAQl, 



(2) 



(3) 



[dQrdQL 

and the sum in (j2]l runs over a two-dimensional grid of < 
15 MeV/c and \ Qj\ < 5 MeV/c, with bin centres located at val- 
ues (Q'j, Q^j) and uniform bin size AQt - AQl - 0.5 MeV/c. 
The Fa and Fb functions describe the A2„ signal and the A^c + 
A^nC + A^acc three-fold background, respectively; cr^ and ctb are 
their statistical errors. The analysis is based on the parametriza- 
tion of Fa and Fg and the precise Monte Carlo simulation of the 
detector response. 

The Fa signal has been simulated fTT] [18] according to an 
accurate model of A2„ production, propagation lfT4l . and inter- 
action with the target medium llT5l[T9H2T| . 

In the background F^, the A^nC and the N^^^ double differ- 
ential spectra were parametrized according to two-body phase 
space and Lorentz boosted to the laboratory frame using the 
observed pion pair spectra IfTTI . The spectrum of A^c pairs is 
enhanced at low Q with Q defined at the point of production, by 
the Coulomb interaction according to the Gamow-Sommerfeld 
factor 

27TM„a/Q 



Ac(Q) = 



(4) 



1 -exp(-27rM;ra/0 
The finite size of the production source and final-state interac- 
tion effects have been calculated ll22l |23]| and applied to sim- 
ulated atomic and Coulomb pairs. An additional momentum- 
dependent correction has been applied to the simulated A^c 
spectrum to take into account a small (< 0.5%) contamination, 
measured by time-of-flight 1*24^, due to misidentified K^K~ 
pairs. Small admixtures of misidentified pp and residual con- 
tamination from e^e^ pairs have been measured and produce 
no effect on the final result. 

The fraction of accidental pairs in Fg was measured by 
time-of-flight to be Wacc - 12.5%, averaged over the pair mo- 
mentum and the different data sets. 
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Figure 1: fit projections of the ff^w" spectrum from data (dots) and sim- 
ulation (MC lines). The top plot shows the experimental spectrum compared 
with the simulated background components (no pionium signal), with (solid 
line) and without (dotted line) Coulomb pairs (At)- The bottom plot shows 
the experimental \ Ql\ spectrum after background subtraction and the simulated 
pionium spectrum. 
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Figure 2: Qr fit projections of the rr^n^ spectrum from data (dots) and sim- 
ulation (line). The left plots show the comparison between the experimental 
spectra and the full simulated background. The plots correspond to different 
Ql regions: top left plot in the signal region (low \Ql\) and bottom left 
plot away from it (higher \Ql\)- The right plot shows the Qt spectrum after 
background subtraction and the simulated pionium spectrum. 



The experimental resolutions on the momentum and open- 
ing angle must be accurately simulated in order to extract the 
narrow pionium signal. Multiple-scattering in the target and 
the spectrometer is the primary source of uncertainty on the 
Qt measurement. In order to achieve the desired Qt reso- 
lution, the scattering angle must be known with ~ 1% preci- 
sion, which is beyond the currently available GEANT descrip- 
tion ||25l . An improved multiple-scattering description was im- 
plemented based on dedicated measurements of the average 
scattering angle off material samples |26|. A cross-check with 
the standard GEANT description was made by comparing the 
momentum evolution of the measured distance between tt^ and 
TT^ at the target |27|. 

The Ql resolution was checked using A decays with small 
opening angle. The widths of reconstructed real and simulated 
A — > pn^ were compared. A 3.4% relative difference was ob- 
served and attributed to residual fringing magnetic field eff'ects, 
multiple scattering in the downstream vacuum channel exit win- 
dow, and to a small misalignment between the spectrometer 
arms. Such effects have been altogether absorbed into an addi- 
tional Gaussian smearing term, of width 0.66 ■ 10"^, convoluted 
with the simulated momentum resolution function. 

The only free parameters in (|2]i are the number of detected 




Figure 3: Coulomb subtracted two-pion correlation function measured in the 
(Ql, Ql) plane, showing the pionium signal, gx is the signed projection of Q 
into a generic transverse axis (azimuthal invariance is ensured by the absence 
of beam and target polarization) 
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atomic pairs (nj"^) and the fraction of non-Coulomb/Coulomb 
pairs (N^^q/N^q^)- The minimization is performed in two- 
dimensional space \Ql\ < 15 MeV/c, Qt < 5 MeV/c, for values 
of the total pair momentum p between 2.6 and 6.8 GeV/c 1281 . 
A constraint on the total number of reconstructed prompt pairs 
is applied such that A^p,.(l - w^cc) - N'c" + ^nC + "a''- 

In Figs. [T] and |2j the \Qi\ and Qt projections of the experi- 
mental prompt n'^Tt' spectrum are shown in comparison to the 
fitted simulated background spectrum (Fa - 0). After subtrac- 
tion of the Fb background, the experimental signal emerges 
at small values of \Ql\ (Fig. [T]i and Qt (Fig. [2| and can be 
compared with the simulated F^ signal. As expected, multiple- 
scattering in the target and upstream detectors broadens the Qt 
signal shape. This is clearly shown in the 2-dimensional plot of 
Fig. [3] The overall agreement between the best-fit experimen- 
tal and simulated spectra is excellent, over the entire Qt, Ql 
domain. 
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Figure 4: The dependence of the measured Ftp, averaged over all data sets, 
from the pionium laboratory momentum and the Monte Carlo prediction cor- 
responding to the ground-state lifetime of 3.15 ■ 10^'^s obtained from the best 
fit. 



4. Pionium breakup probability 

The pionium breakup probability, f br, is defined as the ra- 
tio ha/Na between the number «a of observed pairs from pi- 
onium ionization caused by target atoms and the total number 
A^A of pionium atoms formed by final-state interaction. The lat- 
ter can be inferred by quantum mechanics from the number of 
Coulomb-interacting pairs measured at low Q according to the 
expression ifTOl 



7i:*(Q), 



(5) 



where Q is the domain of integration \Qi\ < 2 MeV/c and 
Qt < 5 MeV/c, yielding ^T* = 0.1301. Differences in detector 
acceptance and reconstruction efficiency for n\ and A^c pairs, 
fA and ec respectively, are taken into account by correcting the 
theoretical factor /T* as 



eA(n) 



(6) 



Those differences arise mainly from the lesser resolution of the 
upstream detectors for identifying close tracks at very low Qt- 
This occurs more frequently for atomic pairs than for Coulomb 
pairs. 

The breakup probability is thus determined as 



«A 
A^A 



1 



A?if<=(Q) 7i:<=='P(Q) 



(7) 



The momentum-dependent /T'^'^p factor (|6]) has been cal- 
culated from fully reconstructed Monte Carlo atomic and 
Coulomb pairs. Using (j6]) and (j?]), 35 independent Pbr val- 
ues are obtained for the five independent data sets and for 
seven 600 MeV/c wide bins of the A2jt momentum from 2.6 to 
6.8 GeV/c, by appropriately folding the momentum dependence 
of 

In Table [T| the fitted yields are given for the different 
momentum-averaged data sets. Overall, more than 2 ■ 10"^ 
atomic pairs have been detected. The reported f br values are 



only indicative of the amount of variation expected with respect 
to the different experimental conditions, and they are not used 
in the final momentum-dependent fit. 

A slight increase of the measured Pbr with increasing pio- 
nium momentum is observed in Fig. |4] (data points), which is a 
consequence of the longer decay path, and hence the greater 
breakup yield, expected at higher atom momenta. The con- 
tinuous curve represents the predicted evolution of Pbr with 
pionium laboratory momentum, for the value of the pionium 
ground-state lifetime t = 3.15 ■ 10"'^ s obtained from this anal- 
ysis. 

The dependence of the breakup probability on the spe- 
cific choice of the integration domain Q has been verified. The 
measured Pbi , averaged over the data sets, is indeed very stable 
versus variations of the \ Ql\, Qt integration limits as shown in 
Fig.|5] 
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Figure 5: Stability of the average P\y, with respect to variation of the: (top) \ Q[} 
(for Qt <5 MeV/c ) and (bottom) Qt (for \Ql\ < 2 MeV/c) integration Umits, 
in 0.5 MeV/c bins. 
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Table 1: Fit results for Qj < 5MeV/c and 16^1 < 15MeV/c. 



Ni, Pbeam 


;r'/ndf 


"A 






^ [ICC 


nr 


94 /zm, 24 GeV/c 
98 //m, 24 GeV/c 
98 //m, 20 GeV/c 
combined samples 


2127/2079 
4288/4149 
4257/4144 


6020±216 
9321+274 
5886±210 
21227±407 


546003+4549 
828554+5811 
496820+4441 
1871377+8613 


45624+4501 
93148+5754 
60867+4397 
199639+8526 


63212+208 
98499+255 
59392+144 
221103+359 


0.441+0.018 
0.452+0.015 
0.472+0.020 



5. Results and systematic errors 

A detailed assessment of the systematic errors affecting the 
Pbr measurement has been carried out, considering all known 
sources of uncertainty in the simulation and in the theoretical 
calculations. The largest systematic error comes from a ~ 1% 
uncertainty in the multiple-scattering angle inside the Ni tar- 
get foil which induces a +0.0077 error on f br. The momentum 
smearing correction can increase Pbr by ~ 2% and thus pro- 
duce a +0.0026 systematic error The double-track resolution 
at small angles can change fbr by 1.1% and generate a sys- 
tematic error of ±0.0014. The admixture of K^K^ changes Pbr 
by ~ 1%. The uncertainty on such contamination is 15% and 
produces a systematic error of +0.0011 on Pbr. The finite-size 
correction to the point-like approximation creates a maximum 
0.8% variation of the simulated yield of Coulomb pairs and a 
systematic error of +0.0011 on fbi. The influence of the final- 
state strong interaction on the r dependence of Pbr is negligi- 
ble ifTSl I^?!. The trigger response efficiency was measured us- 
ing minimum-bias events and accidental pairs from calibration 
runs. The efficiency is high and quite uniform in the selected 
Qt, Ql domain and it drops by ~ 2% per MeV/c at \Qi\ > 15 
MeV/c. The simulated and experimental trigger efficiencies 
agree to better than 0.5%, in the same \Qi\ range. This max- 
imum deviation increases the breakup probability by ~ 3% and 
thus produces a systematic error of +0.0004. Background hits 
in the upstream spectrometer region, generated by beam and 
secondary interactions in the target region, are the source of a 
+0.0001 systematic error on Pbr- The effect of the lower purity 
of the 94 ixva Ni target foil compared to the 98 jim is an under- 
estimation of fbr by ~ 1.1%. This corresponds to a systematic 
error of +0.0013 for the corresponding data set. 

The dependence of Pbr on the atom lifetime t, its momen- 
tum, and the target parameters has been extensively studied for 
several target materials, both by exactly solving the system of 
transport equations lfT4l [TSl describing the excitation/de- 
excitation, breakup and annihilation, and by simulating fT5\ 
the A2n propagation in the target foil. The precision reached 
by these calculations is at the level of 1% |291, which is re- 
flected in a +0.0042 systematic error on f br for a lifetime r = 
3.15 ■ 10"'^ s. The result of these calculations defines three 
functions Ph^ij, p), one for each of the combinations of target 
thickness and beam momentum. The functions P\,r(T,p) are 
further convoluted with the experimental momentum spectra 
of Coulomb pairs inside the seven (600 MeV/c wide) mom- 
entum slices of the pionium laboratory momentum, from 2.6 
to 6.8 GeV/c. This approach ensures that within each slice the 
non-linear dependence of Pbr(T") on the laboratory momentum 
is negUgible. 



Coulomb pairs, which have a momentum spectrum similar 
to that of atomic pairs, are taken from prompt pairs in the Q 
region away from the A2„ signal, after subtraction of the non- 
Coulomb contribution. The values of the systematic errors are 
summarized in Table |2l 



Table 2: Summary of systematic errors on Ftp. 



source 


(T 


multiple scattering 


+0.0077 


momentum smearing 


+0.0026 


double-track resolution 


+0.0014 


K^K- and pp 


+0.0011 


trigger simulation 


+0.0004 


background hits 


+0.0001 


target impurity 


+0.0013 


finite size 


+0.0011 


calculation of PhAT) 


+0.0042 


Overall error 


+0.0094 



6. Conclusions 

Finally, the Pbr measurements, obtained for the different ex- 
perimental conditions and A2„ momentum ranges, and their pre- 
dicted P\,i-(t, p) values (see Fig|6]), were used in a maximum 
likelihood fit of the lifetime t [30|. Both statistical and system- 
atic uncertainties were taken into account in the maximization 
procedure. 




X, s 



Figure 6: Function Pbi ir) corresponding to the dependence on pionium lifetime 
of the breakup probability for different targets. 
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Our final measurement of the ground-state A2;r lifetime yi- 

eldsr.(3.15tf>Lt?^Lj><10-'^s. 

Taking into account A2jt — > yy and using formula ([TJ, we 
obtain the nn scattering length difference 

U „ I _ /n oc^^+O.OOSOl +0.00781 \ .^-1 /ox 

|ao - ail - ^0.2533_oo,)78|^j^^ -o.oovBlsystj ' i^) 

where the systematic error includes the 0.6% uncertainty in- 
duced by the theoretical uncertainty on the correction 6. 

In conclusion, we have measured the ground-state lifetime 
of pionium with a total uncertainty of ~ 9%. This represents the 
most accurate lifetime measurement ever obtained and has al- 
lowed us to determine the scattering length difference |ao - ail 
with a ~ 4% accuracy. Our result is in agreement with val- 
ues of the scattering lengths obtained from Ke4 |8J and K3„ f9l 
decay measurements using a completely different experimental 
approach. 
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